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ABSTRACT 

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact 
sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due 
to confusion in a background of fluctuating surface brightness." The histogram of the flux densities of extracted sources shows 
a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This 
empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus 
noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with 
statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used 
extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric 
operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 
to 500 yum within Galactic plane fields at / = 30° and / = 59°. We find that the exponent of the power spectrum is about -3. At 250 fim, 
the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus 
noise scales linearly with the median brightness. For a given region, the wavelength dependence of the amplitude can be described 
by the square of the spectral energy distribution (SED) of the dust emission. Generally, the confusion noise will be a worse problem 
at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower 
spatial frequencies, but the photometric signal to noise will also depend on the relative SED of the source compared to the cirrus. 

Key words. ISM: general, structure — Stars: formation, protostars — Submillimeter: ISM 



""j 1. Introduction the Galactic plane, which is also bright. Cirrus noise varies 
^ with cirrus brightness (Sect. |3]l, introducing further complex- 
H ; Cirrus noise, which is operationally defined as the "statisti- ity to the problem. Cirrus fluctuations characteristically decrease 
>^ , cal error to be expected in photometric measurements due to ^jth decreasing spatial scale (Sect. [3]l, but even with the im- 
confusion in a background of fluctuating surface brightness" pj-oygd angular resolution of Herschel, cirrus noise remains a 
( |Gautieretal. L 11992D, is a major issue limiting the cataloging of dominant factor. The Herschel observation planning tool HSpot 
compact sources that underpins the study of the early stages of (www.ipac.caltech.edu/Herschel/hspot.shtml) has a built-in con- 
star formation in the interstellar medium. Examination of wide f^gjon ^^^^^ estimator to provide "on-line guidance on where to 
range of mass of the stellar precursors, requires measurement g^pect fundamental detection limits for point sources that can- 
of sources with a wide range of luminosity, or at each wave- ^^^ ^e improved by increasing the integration time." With such 
length, flux density. Stars form where there is abundant mate- guidance the serious impact of cirrus noise was anticipated, and 
rial, and so the ciiTus brightness in the field tends to be high, ^s shown below it can now be quantified directly using submil- 
Furthermore, many studies, both targeted and unbiased, are in limeter data 



* Herschel is an ESA space observatory with science instruments To appreciate the problem at its fundamental level con- 

provided by European-led Principal Investigator consortia and with im- sider actual catalogs of extracted sources (Mol inaril 120101: 
portant participation from NASA. iMolinari et all 1201 Oah in two degree-sized Hi-GAL Galactic 



p. G. Martin et al.: Direct Estimate of Cirms Noise in Herschel Hi-GAL Images 



20 ■ 



5 lol- 

E 



_ _l ' 






^- 



Flux density (Jy) 

Fig. 1. Solid histogram: 136 sources cataloged at 250 fim in a 
representative large degree-sized sub-region of Hi-GAL / - 30° 
field (f30: see Sect.|2]). Dashed histogram: 70 sources cataloged 
in a sub-region of Hi-GAL / = 59° field (f59) that has a factor 
5.8 lower median brightness. 



plane fields jMolinari et all 1201 Obh (see Sect.|2]l. The solid his- 
togram in Figure [U shows the flux densities of sources cataloged 
at 250 yum in the brighter field. Note the falloff in source counts 
at flux densities less than 10 Jy. This falloff is not an intrinsic 
property of the underlying population of sources. We show be- 
low that this falloff is the expected consequence of cirrus noise, 
which can be quantified independently of the making of the cat- 
alogs. In the fainter field (dashed histogram), the cutoff is lower 
because of decreased cirrus noise. Clearly, it will be important to 
account for the varying dramatic effects of cirrus noise in order 
to recover the statistics of the intrinsic faint-source population. 

Cirrus noise is operationally defined and so for a partic- 
ular measuring strategy, such as fitting compact sources with 
Gaussians as used in Hi-GAL, this can actually be estimated di- 
rectly from the source-removed maps (Sect. [3]). There is also an 
extensive literature on "estimating confusion noise due to ex- 
tended structures giv en so me statistical properties of the sky" 
(Gautier et al.', '1992'); see (Kiss et al.', '2003'; Jeong et all l2005t 
[Miville-Deschenes et al., 2007; Roy et al., 2010). The cirrus 
brightness statistics are well described by a power spectru m, ap- 
propriate for Gaussian random fields. Gautier et al.' (1 19921) found 
some non-Gaussianity and [Miville-Deschenes e t alJ (l2007h 
found non-vanishing skewness and excess kurtosis in the under- 
lying brightness fluctuation fields. Nevertheless, these are not 
large effects and for estimating the variance the power spectrum 
is demonstrably still a powerful tool, particularly because of the 
insight it provides into how the cirrus noise depends on vari- 
ous statistical and operational parameters. Therefore, we con- 
nect the new direct estimates with what can be obtained using 
power spectra (Sect. [3]). In Sect. |4] we show that reliable power 
spectra can be obtained even from first-generation images pro- 
cessed for Hi-GAL. Finally, we return to how estimates of cirrus 
noise should ultimately explain the faint-end cutoff in forthcom- 
ing source catalogs (Sect.|5]). 

2. Observations 

Two fields were observed (at / ~ 30° an d / ~ 59°) during SDP as 
part of Hi-GAL ( Mohnari et all 12010^ . using the Paraflel mode 
with both nominal and orthogonal scanning to acquire data in the 
PACS 70 and 160//m bands and all three SPIRE bands (250, 350 
and 500 jum). The data processing using HIPE and map-making 



using ROMAGAL are described by iMolinaril (1201 Ol) and refer- 
ences therein. The maps were converted into brightness units 
I MJy/ sr) and the DC offsets recovered as described by Bernarg 
20101) . The offsets are not needed for the cirrus noise analysis, 
but are needed to describe, for example, the median brightness 
of the map used in Sect. [3] and Sect.|4] 

Compact sources have been detected and then quantified us- 
ing a Gaussian model (Molinari, 2010; Molinari et al., 2010af). 
Using the modeled properties, the sources have been removed 
to produce the images used for the analysis here. There remain 
some bright, fairly compact peaks not meeting the classification 
criteria for compact sources, but probably gravitationally influ- 
enced and not standard cirrus structure. There are also some arte- 
facts from removal of the brightest compact sources. Given the 
preliminary nature of these source-subtracted images, for this 
initial analysis we have analysed large degree-sized sub-regions 
in which these effects are minimized. They are otherwise repre- 
sentative of bright cirrus. The brightest is from the fairly homo- 
geneous low-longitude half of the / ~ 30° map, delineated by 
29 < I < 29.83° and -0.5 < b ^ 0.33°, with median bright- 
ness < /250 > = 1390 MJy sr '. For convenience, we refer to 
this field as f30. A contrasting region, in having a lower median 
brightness 241 MJy sr"', is the lower right quarter / ~ 59° map, 
delineated by 58 ^ Z ^ 59° and -1 < Z? < 0° (denoted f59). 

3. Direct estimates of cirrus noise and relationsliip 
to power spectra 

In the literature, cirrus noise is defined operationally for a pho- 
tometric measurement template used to evaluate the flux density 
of a compact source on a fluctuating background (not simply 
"finding" it). We obtain a direct estimate of cirrus noise numeri- 
cally by placing the template randomly in a cirrus map and find- 
ing the rms of the apparent " source" flux dens ity, o-jii-i-us. In the 
context of power-law cirrus. iGautier et all d 1992 ) quantified this 
analytically for several measurement templates, and these plus 
simulated cirrus maps have been used to validate our numeri- 
cal approach in detail (see below). Among these templates is the 
"apert ure plus reference annulus" c hosen by Helou & BeichmanI 
(119901) and adopted by Kiss et alJ (1200 ih and lRovetal] (1201(1 ? 
In practice, compact sources are often fit with a model template 
consisting of a Gaussian and a planar inclined plateau, with a 
footprint larger than the extent of the Gaussian to provide a ref- 
erence area to characterize the "background." We adopt a radius 
of 1.82 times the FWHM of the Gaussian used. This Gaussian- 
based strategy has the practical advantage of still being useful 
for moderately crowded sources, unlike simple aperture plus 
annulus photometry. As a specific and relevant example using 
the 250jum Hi-GAL fields, the source FWHM is typically 1.5 
times larger than the diffraction-limited PSF and so we adopt a 
Gaussian template with this larger FWHM. We also assume that 
the source-subtracted map is a good proxy for the cirrus in the 
field. The estimated cirrus noise is 1.7 Jy for f30 and 0.23 Jy for 
the fainter f59, scaling roughly as the median brightness which 
changes by a factor 5.8. 

As mentioned, there is extensive literature in which the cir- 
rus noise is quantified using the power spectrum. Furthermore, 
these authors showed that the power spectrum of Galactic cir- 
rus follows a power law Pdrmsik) - Piko)(k/ko)", quantified by 
an amplitude Pq = P( ko) at some fiducial kg, and an exponent 
a that is typically -3 (Miville-Deschenes et all l2007h . as also 
found here. We adopt ko = 1.0 arcmin"', a scale close to the 
beam sizes, to avoid issues of extrapolation if a is not quite -3. 
The essentials of the detailed JGautier et alJ(ll992h analysis of the 
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"cirrus noise" for a telescope with mirror diameter D working at 
wavelength A and a - -T> can be summarized as: 



Table 1. Parameters" for SPIRE T{k) at 250, 350, 500^^^ 



O'ciiTus — ^t "f 



,2.5 



2.5 



Pq 



( A/250 fim 

\ D/3.5 m / \lO''Jy 



sr 



0.5 



Jy. 



(1) 



Here Rt is a dimension of the measurement template (e.g., aper- 
ture diameter or FWHM of Gaussian) in units of A/D (explicitly 
in the formula with the same power). If the compact sources be- 
ing measured are extended, as will be the case in many Galactic 
surveys including those illustrated here, then R, needs to be cor- 
respondingly larger. The consequent increase in o-dn-us is cap- 
tured by the factor Rf^ for the range of interest. A, is an am- 
plitude appropriate to a particular measurement template. A, ac- 
tually changes significantly with R, for small R,, but this is not 
relevant because no practical application would use photomet- 
ric templates smaller than the diffraction limit; for the range of 
in terest, there is only a weak /?P dependence (see also Fig. 3 
in lGautier et al.L[l992l) . Neither A, or R, depend on wavelength; 
that dependence is carried explicitly in the A/D term and implic- 
itly in the square root dependence on Pq. 

We have made simulated power-law cirrus maps with a - -3 
and directly evaluated the cirrus noise with various templates 
for a range of/?,. We recover the predicted /?p dependence and 
specific values of A,, like 0. 034 predicted for the wid ely-used 
aperture-annulus template of lHelou & BeichmanI (1 19901) with R, 
near 1.6. From these simulated maps we also derived the ampli- 
tude appropriate to the Gaussian fitting template for the empir- 
ical range of Rf near 1.8 for the actual sources. When the sim- 
ulated cirrus map has resulted from convolved with the actual 
Neptune beam (Sect.©, A, - 0.054 (A, - 0.065 for the nominal 
Gaussian PSF). 

4. Properties of observed power spectra 

Details of computing t he power spectrum and its errors may be 
found, for example, in iRov et al.l (l2010l) . In practice, contribu- 
tions to the total power spectrum come not only from diffuse 
dust emission, but also from point sources, the cosmic infrared 
background (CIB), and "noise" (which might reasonably be as- 
sumed to be fairly white). When these components are statisti- 
caUy uncorrelated, the total pow er spectrum can be expressed as 
dMiville-Deschenes et ani2007h : 



P(k) = r(^) [Feirrus(^) + ^source(^) + Pcisik)] + N(k). 



(2) 



For bright Galactic plane fields, the contribution from the CIB to 
the power spectrum is negligible. We have removed the brighter 
compact sources, but others at or below the detection threshold 
must remain and need to be accounted for by Psource- While the 
"noise" A^ is measurable, and of interest, it too makes an insignif- 
icant contribution. 

r{k), the power spectrum of the PSF, de- 
cays at large k. For the SPIRE bands we 
used the empirical PSF from scans of Neptune 
( |ftp://ftp.sc iops.esa.int/pub/hsc-calibration/SPIRE/PH OT/Beams ) 
and co mputed r(k) using the technique developed bv lRov et al.l 
( I2OIOI) . A Gaussian approximation is poor and adversely 
impacts the extraction of Pcimis(fc), but the fitted functions in 
Table [T] are good to within a few percent for the range of interest 
k < k,„ax/2, where k„ax is the highest spatial frequency in the 
power spectrum for these approximately Nyquist-sampled maps 
(5 arcmin"' at 250;um). 

The top set of curves in Figure |2] shows the power spectrum 
at 250 //m for the f30 field. Compared to the power spectra from 
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Fig. 2. Top curves: power spectrum at 250 jum for the f30 field 
(Sect. |2]i. Decay at large k in lower of the two curves, P(k), 
is due to r(k) from the PSF. Plateau at high k is from resid- 
ual "noise" A^ in the map, 2.4 x 10'* Jy^ sr"'. At very small 
scales, where the astronomical signals become coiTelated within 
the beam, P meets A^. Upper curve is Pcinus(^)+Psoui-ce, after sub- 
tracting the noise and dividing by r(k) to remove the effect of the 
beam. Upper curve has been fit over the range 0.06 arcmin ' < 
k < 2.5 arcmin"' by a simple model (dotted curve) consisting of 
a power law (dashed line) with Pq - (58+1) x 10^ Jy^ sr"' 
plus a constant Psomce = (4.9 + 0.3) x lO'' Jy^ sr"' (dash dot 
line). The power-law exponent is -2.74 + 0.03. Bottom curves: 
power spectrum for the fainter f59 field. The "noise" level, Po, 
and Psource have all decreased, roughly as the square of the me- 
dian brightness, to 7.3 x 10^ and (1.13 + 0.03) x 10* and 
(8.6 + 0.6) X 10* Jy^ sr"', respectively. The power-law expo- 
nent -2.81 + 0.03 is not significantly different. 



BLAST shown bv lRovetaJ] (1201 Ol) . Pik) can be determined to 
larger k, as expected because of the smaller Herschel beam. The 
lower set of curves shows similar results for the fainter f59 field. 
The steps in fitting the data by equation |2] to find PciiTus(^), and 
the results of the fit, are described in the figure caption. The fit 
is very good (reduced ;t-^ ~ 1.5). 

For the fainter f59 field, the noise level, 7.3 x 10^ Jy^ sr"', 
is comparable to the detector noise predictable by HSpot. It is 
in fact slightly lower than for the much fai nter Polaris field ob- 
serve d with the same coverage (see Fig. 3 in lMiville-Deschenesl 
I2OIO ). Tests indicate that this slight improvement comes from 
using ROMAGAL rather than the naive map-maker in HIPE. 
However, note that the noise level for the brighter sub-region 
is dramatically higher, 2.4 x 10"* Jy^ sr"'. This appears to scale 
as the square of the median brightness of the map and so must 
measure residual map artefacts, not detector noise. Nevertheless, 
these are already excellent images with the "noise" small com- 
pared to Pcimis(^) over the range of interest k < 2.5 arcmin"'. 
The presence of Psource is only subtly evident in a slight con- 
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Fig. 3. SED obtained from the square root of the ampHtudes Pq 
obtained from fits to power spectra in the f30 field. SoHd curve 
is from the fit of a yS = 2 modified blackbody, with temperature 
23.6 K. 



cave curvature of the power spectrum at large k and is adequately 
modeled there as a constant. 

For f30 we find a = -2.74 + 0.03 and in the fainter f59 
field -2.81 + 0.03, not significantly different. These are close 
to estimates in the Galactic plane bv lRov et al.l (1201 Ol) . and also 
remarkably similar to what has been found at 100/zot for high- 
latitude diffuse cirrus jMiville-Deschenes et al.L l2007h . despite 
quite different conditions and geometries which could in princi- 
ple aff'ect the cirrus fluctuations. This supports the adoption of 
a = -3 for equation[T]to quantify CTcimis in the submillimeter. 

At 250 /im we have examined two further regions, the up- 
per right quarter of the / - 59° map and the Polaris field. 
From this we conclude that Pq scales roughly as the square of 
the median brightness. A corollary is that a-^-i^^^^ should scale 
roughly as the median brightness, as we have already found di- 
rectly (Sect. [3]). This behavior is as though the inherent statistical 
structure in brightness ly is the same, just scaled up and down. 
Empiric ally, at lOOyum a steeper scaling has been seen for bright 
cirrus (.Miville-Deschenes et all 12007.) . and so this should be re- 
visited for all wavelengths when more fields are available. One 
factor that is changing from field to field, and within fields, is the 
dust temperature T^ (Bernard, 2010). Even with the same dust 
column density structure, this will modulate the cirrus bright- 
ness, with stronger effects at 100 jum than in the submillimeter 

For a given region, images of optically thin dust emission 
at different passbands could be simply related by the scale fac- 
tor KyBy{Ti) or the (relative) spectral energy distribution (SED). 
Here, Ky is the dust emissivity. Because of averaging along the 
line of sight and over different grain components with potentially 
differe nt T^, this is a simplification. Nevertheless, iRov et all 
(12010') demonstrated with BLAST and IRAS data how an SED 
with a reasonable T^ can be recovered from the wavelength de- 
pendence of Pj (conversely, Pq varies as SED^). 

We have examined this with the Herschel data, measuring 
the power spectra and finding Pq for each of the five passbands. 
The results for the / ~ 30° sub-region are plotted in Figure [3] 
These data were then fit with a single-temperature modified 
blackbody {Ky oc y^ or/? = 2). We find that the functional form 
is satisfactory, and that T^ = 23.6 + 1.0 K. The temperature has 
been fit independently pixel by pixel (Bernard, 2010). For this 
sub-region, with/? - 2 the average temperature is 22.5 K (range 
16 - 28.5 K, median: 22.5 K). The close agreement is reassur- 



ing. Thus, other things being equal, the SED of CTcinus will be 
close to that of the cirrus emission. In measuring the flux density 
of a compact source, the wavelength dependence of the signal 
(source SED) to noise (diffuse dust SED) will depend on the rel- 
ative temperatures as well. Sources cooler (hotter) than the cirrus 
will be detected at a relatively higher (lower) S/N at long wave- 
lengths. However, equation [1] also shows that in practice cirrus 
noise will be more severe for the long-wavelength bands because 
of the larger beam, offsetting this effect on S/N for cooler sources 
and compounding it for hotter ones. 



5. Cirrus noise and faint-end cutoffs in catalogs 

We can estimate the cirrus noise from the power spectra and 
equation!!] For f30 and f59 at Pq is 58 and 1.1 x 10^ Jy^ sr', 
respectively. When fit by Gaussians, sources in these Hi-GAL 
fields require typically R, = 1.8. A, = 0.054 has been evaluated 
from simulations for this range of/?,. We find Cronus = 1-8 and 
0.25 Jy, respectively. These compare favorably with the values 
estimated directly, 1.7 and 0.23 Jy, respectively (Sect.[3]l. In the 
f30 field, a 5-cr catalog would be cut off at about 9 Jy. The fainter 
f59 field is also affected, but would have a lower 5-cr cutoff of 
1.3 Jy. Cirrus noise therefore severely limits the depth of the cat- 
alogs and generally in Hi-GAL Galactic plane fields will be the 
pre-dominant factor When catalogs from multiple fields with 
different median brightnesses are combined/interpreted, account 
must be made for the differing impact of cirrus noise. 

The faint-end cutoffs in the catalogs of compact sources il- 
lustrated in Figure[T|are close to the estimated influence of cirrus 
noise. Quantitat ively, th is must be somew hat of a coinci dence. 
As described by MoUnari (2010) and MoH nari et"an ('2010a). the 
catalog is not selected on the basis of signal to noise (S/N) of the 
flux densities. No errors have been tabulated for the flux densi- 
ties from the Gaussian fits. Judging from their a posteriori S/N 
estimator, which involves peak flux (Jy/beam), there is a wide 
range, including some below 5. Simulations to assess catalog 
completeness find that the peak flux for 80% completeness is 
4.1 Jy/beam over the entire / ~ 30° map, and so for the typi- 
cal source FWHM, the corresponding flux density completeness 
limit is 9.2 Jy close to our estimated 5-cr threshold and the ob- 
served histogram peak in Figure[T] For the fainter I ~ 59° survey 
there is a lower 80% completeness limit (1.6 Jy); peak in the his- 
togram occurs at a lower value, again close to this limit and our 
estimated 5-cr threshold. 

At longer wavelengths, cirrus noise will be more limiting for 
both source detection and source flux density determinations, be- 
cause of the A^-^ (beam) dependence. This highlights another im- 
portant consideration for band-merged catalogs, that the S/N for 
cataloged sources will be wavelength dependent. 
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